(Original source: this answer in Mathstackexchange
In $\mathbb{R}^m$ we have a natural basis, the canonical one, and we have a natural volume form, the classical determinant. Indeed, we can think that the choice of this basis gives rise to an inner product (by establishing this basis as orthonormal) whose associated volume form is the determinant. That is, the determinant measures the same volume on parallelepipeds that we would obtain by measuring the edges and angles with the inner product. The volume form is not determined, therefore, by the basis but by a class of basis with the same associated inner product (one basis can be obtained from the other by an orthogonal transformation).
But we want also to compute "intermediate volumes": lengths, areas, 3-volumes,... up to the $m$-volume, and this can be done by using the inner product (classical scalar product). But there is a formula to compute them directly. If $A$ is the matrix whose columns are the $n$ vectors $v_1,\ldots,v_n \in \mathbb{R}^m$ whose $n$-volume want to be computed then the $n$-volume is
$$ D(A):=\sqrt{\det(A^T A)}. $$I think it is called something as Gram matrix.
Observe that this formula includes as particular cases the 1-volume or "length", and the $m$-volume (determinant of a square matrix).
In order to prove this formula, we can consider that the $n$ vectors $v_1,\ldots,v_n \in \mathbb{R}^m$ are linearly independent (if not, their $n$-volume would be zero) and look for an orthogonal transformation $T$ of $\mathbb{R}^m$ (with matrix $U$ in the canonical basis) that sends the $n$-dimensional subspace spanned by them to $span(e_1,\ldots,e_n)$ being $\{e_i\}$ the canonical basis of $\mathbb{R}^m$. Since the transformation is orthogonal, the $n$-volume should be the same, that is,
$$ Vol(A)=Vol(UA). $$But this volume, since is "inside" $\mathbb{R}^n\subset \mathbb{R}^m$ can be computed classically as
$$ \det(trunc(UA)) $$where $trunc$ means selecting the first $n$ rows (the others are null).
And finally, observe that
$$ \det(trunc(UA))=\sqrt{\det(trunc(UA)^Ttrunc(UA))}= $$ $$ =\sqrt{\det((UA)^T UA)}=\sqrt{\det(A^T A)} $$since the truncated part is a 0 block matrix and since $U^T U=Id$.
So $Vol(A)=\sqrt{\det(A^T A)}$.
I we had started with a different inner product $\langle-,-\rangle$ (or with another basis, not corresponding to an orthogonal transformation of the canonical one), I think the formula would be nearly the same:
$$ Vol(A)=\sqrt{\det(A^T G A)} $$where $G$ is the matrix of the inner product, that is, $G=(g_{ij})$ with $g_{ij}=\langle e_i, e_j \rangle$.
Well, indeed we must define first what is the "volume determined by the classical scalar product". Let's call it the metric volume $Vol(A)$. I would do it inductively: the metric volume of $v_1,\ldots,v_n$ is the metric volume of $v_1,\ldots,v_{n-1}$ multiplied by the projection of $v_n$ over the orthogonal complement of $v_1,\ldots,v_{n-1}$ in $v_1,\ldots,v_n$.
With this in mind, the result is evident when the given vectors are multiples of the basis vector. If we assume that they are given in the correct order then the matrix $A=\Sigma$ is a diagonal matrix whose entries are the factors and therefore $D(\Sigma)=Vol(\Sigma)$.
On the other hand, observe that if $U$ is an orthogonal transformation and $D(B)=Vol(B)$ then $D(UB)=D(B)=Vol(B)=Vol(UB)$, where the last equality holds because $U$ preserves angles and distances. So if the given vectors are multiples of the result of an orthogonal transformation (reflection-rotation) applied to the basis vectors of $\mathbb{R}^m$ then $D(A)=Vol(A)$, since $A=U\Sigma$, where $U$ is a orthogonal matrix of order $m.$
In the general case, the matrix $A$ can be seen like the matrix of a transformation $T$ sending the canonical basis of $\mathbb{R}^n$ to the vectors $v_1,\ldots,v_n.$
Taking into account the SVD of this matrix
$$ A=U\Sigma V, $$we can think in $T$ like a sequence of:
We only need to show that $D(\Sigma V)=Vol(\Sigma V)$. We can check that
$$ D(\Sigma V)=\sqrt{\det((\Sigma V)^T \Sigma V)}=\sqrt{\det\Sigma^T\Sigma}=D(\Sigma)=\prod \sigma_{ii} $$But on the other hand, we can think of $\Sigma V$ as a transformation of $\mathbb{R}^n$ into the first $n$ components of $\mathbb{R}^m$, simply by truncating the first $n$ rows. And
$$ \det(\mbox{Trunc}(\Sigma V))=\prod \sigma_{ii} $$is nothing but the expansion factor of the volume of parallelepipeds for this linear transformation. So $Vol(\Sigma V)=Vol(\mbox{Trunc}(\Sigma V))=\prod \sigma_{ii}$
since if we restrict to the first $n$ basis vectors of $\mathbb{R}^m$, $\det(trunc(\Sigma V ))$ is the expansion coefficient, and the original parallelepiped in $\mathbb{R}^n$ had volume 1.
________________________________________
________________________________________
________________________________________
Author of the notes: Antonio J. Pan-Collantes
INDEX: